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ABSTRACT 

Three planets with minimum masses less than 10 M ffi orbit the star HD 40307, 
suggesting these planets may be rocky. However, with only radial velocity data, 
it is impossible to determine if these planets are rocky or gaseous. Here we exploit 
various dynamical features of the system in order to assess the physical properties 
of the planets. Observations allow for circular orbits, but a numerical integration 
shows that the eccentricities must be at least 10 -4 . Also, planets b and c are so 
close to the star that tidal effects are significant. If planet b has tidal parameters 
similar to the terrestrial planets in the Solar System and a remnant eccentricity 
larger than 10 -3 , then, going back in time, the system would have been unstable 
within the lifetime of the star (which we estimate to be 6.1±1.6 Gyr). Moreover, 
if the eccentricities are that large and the inner planet is rocky, then its tidal 
heating may be an order of magnitude greater than extremely volcanic Io, on 
a per unit surface area basis. If planet b is not terrestrial, e.g. Neptune-like, 
these physical constraints would not apply. This analysis suggests the planets 
are not terrestrial-like, and are more like our giant planets. In either case, we 
find that the planets probably formed at larger radii and migrated early-on (via 
disk interactions) into their current orbits. This study demonstrates how the 
orbital and dynamical properties of exoplanet systems may be used to constrain 
the planets' physical properties. 
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1. Introduction 

Mayor et al. (2008) recently announced the discovery of a system of three planets, b, c, 
and d, orbiting the K dwarf HD 40307. This system is unique because it is the first detected 
system in which all three companions have minimum masses less than 10 M e . Moreover, 
the innermost planet, b at 4.2 M e , has the lowest minimum mass yet detected by radial 
velocity methods. In addition to relatively low masses, the system is striking in that the 
planets appear to lie close to a Laplace-like resonance: very small eccentricities ( <0.01) and 
period ratios near 4:2:1. However, as noted by Mayor et al. (2008), the observations rule out 
such a resonance chain with high confidence (> 10o"). 

If, as is most probable, the actual masses are similar to the minimum value, it is natural 
to wonder whether these bodies are larger versions of the rocky planets in the Solar System 
("super- Earths"), or smaller versions of our gaseous planets ("mini-Neptunes"). The only 
currently available method to make a direct assessment of these two possibilities requires 
transit data, but none has been reported, so we must rely on indirect means. One possibility 
is to consider the theoretical modeling of terrestrial and gaseous planet formation. However, 
that approach leads to uncertainty. For example, core-accretion models predict that a solid 
core requires 2-10 M e in order to accrete a hydrogen envelope (Pollack et al. 1996; Boden- 
heimer, Hubickyj, & Lissauer 2000; Ikoma et al. 2001; Hubickyj, Bodenheimer & Lissauer 
2005). Therefore, planet formation models cannot yet constrain the physical nature of the 
4 — 10 M e planets in this system. 

Here we exploit another method for constraining the properties of exoplanets: The 
orbital history since formation. In our Solar System, the rocky and gaseous planets' responses 
to tides are very different; solid, rocky bodies dissipate tidal energy more effectively (smaller 
Q' values) than their gaseous counterparts. Tides result in orbital migration at rates that 
can be orders of magnitude different depending on whether the planet is rocky or gaseous 
(Jackson et al. 2008a, 2008b). Also, the effectiveness of tides falls off rapidly with distance 
such that in multiple planet systems with close-in planets, the innermost planet has usually 
experienced tidal migration, while the others have undergone little. By modeling the past 
tidal evolution of a system, the inner planets' properties may be constrained by forbidding 
past events, e.g. mean motion resonance crossing or crossing orbits, that would have led to 
orbits inconsistent with the current system. For example, Barnes et al. (2008) considered the 
GJ 581 system (Udry et al. 2007) and showed that planet c ( >5 M e ) could be terrestrial, 
but cannot have the same tidal parameters as the present-day Earth. Jackson et al. (2008b) 
considered the GJ 876 system (Rivera et al. 2005) and showed that planet d ( >7.5 M e ) 
cannot be terrestrial because, considering tides, e would have been ps 1 less than ~ 30 Myr 
ago, with internal heating rates up to 10 5 times that of Io. 
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For the HD 40307 system, we can exploit the proximity of mean motion resonances to 
constrain the tidal evolution of the innermost planet, b. We focus on two end-member cases 
for HD 40307 b: rocky or gaseous, but allow for other possibilities. To exploit dynamical 
constraints, we also use what we know about the system's age, as well as how it may have 
formed. As we show in the following sections, the possibility that the planets are terrestrial 
in character (i.e. "super-Earths") seems unlikely, but cannot be ruled out. However, the 
gaseous case (i.e. "mini-Neptunes" ) is less constrained. Nor can we rule out an intermediate 
case, that is tidal parameters in between those of the terrestrial and gaseous planets in our 
Solar System. In § 2 we describe our dynamical and tidal models. In § 3 we present our 
results for this system. In § 4 we infer the character of these planets, and identify likely 
formation scenarios. In the Appendix we estimate the age of the star, and, by extension, the 
system. 



2. Methods 
2.1. Planet-Planet Interactions 



We will consider the oscillations of the planets' orbits with the N-body code HNBody[j], 
which includes general relativistic effects. For these integrations, we require numerically- 
induced energy changes to be less than one part in 10 4 , which is adequate precision to 
produce reliable results (Barnes & Quinn 2004). 



2.2. Tidal Evolution Models 

For our tidal model we use conventional equations assembled by Goldreich & Soter (1966; 
see also Jackson et al. 2008a; Barnes et al. 2008; Ferraz-Mello et al. 2008). The evolution of 
semi-major axis a and eccentricity e (to second order in e) can be modeled as 

da _ f 63 s/GM!R^ 2 + 9 y^%\ _ 11/2 (1) 



dt \2 m p Q' p 2 
de _ (63 y/GM*B* 225 y/G/M^,,,,, \ _ |: . ., 



,a- Li/2 e, (2) 
dt V 4 m p Q' p 16 Q; J y ' 

where G is Newton's gravitational constant, M* is the stellar mass, R p is the radius of the 
planet, m p is the planet mass, Q' and Q'^ are the planet's and star's tidal dissipation function 



1 Publicly available at |http: / /www.astro.umd.edu/ ~rauch /HNBody 
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divided by two-thirds its Love number, and R* is the stellar radius. In Eqs. (p]-[2]) the first 
terms represent the tide raised on the planet by the star, and the second terms the tide raised 
on the star by the planet. This model assumes the tidal components maintain a constant 
phase lag from the line connecting the centers of mass of the two bodies, and is consistent 
with observations of the Galilean satellites of our Solar System (Aksnes & Franklin 2001). 
See Jackson et al. (2008a) or Barnes et al. (2008) for more discussion of this model, but note 
that other plausible models, also consistent with observations in the Solar System, exist (see 
e.g. Hut 1981; Neron de Surgy & Laskar 1997; Eggleton et al. 1998; Mardling & Lin 2002; 
Efroimsky & Lainey 2007; Dobrovolskis 2007; Ferraz-Mello et al. 2008). 

Our model does not include the effects of secular interactions between the planets (see 
e.g. Mardling 2007). Such effects could potentially play a role, but we will show that our 
conclusions regarding the HD 40307 b system are probably not impacted by the neglect of 
this effect. If it were included, then the timescale for eccentricity evolution would likely be 
slower, and the tidal evolution of the innermost planet would change the eccentricities of 
other planets in the system. (Exterior planets modify the eccentricity of the inner planet, 
as tidal evolution damps it.) Eqs. ([I]) and (J2J) allow a reasonable description of the tidal 
evolution of the system. 

We may also determine the amount of heat generated in a body due to tidal friction: 

63 (GM.)WM.Bj / 
H= T ^ a e (3) 

(Jackson et al. 2008b). H represents the internal heating rate, but for geophysical consid- 
erations it is useful to express the heat as a surface flux h = H/47T 2 R p . For reference the 
heat flux on the Earth (due to radiogenic processes) is 0.08 W/m 2 (Davies 1999), Io's is 2 
- 3 W/m 2 (Peale et al. 1979; McEwen et al. 2004), and Europa's, scaling from Io's, could be 
~ 0.2 W/m 2 (O'Brien et al. 2002). 

A key parameter is Q' , which parameterizes the planet's tidal response to the star. In 
principle Q' p may take any value, but in our Solar System rocky and gaseous bodies tend 
to cluster around two Q' p values separated by several orders of magnitude. For terrestrial 
bodies Q' p ~ 500 (Dickey et al. 1994; Mardling & Lin 2004; Lambeck 1977; Yoder 1995). 
For gaseous bodies it is common to adopt Q' p ~ 10 5 (Banfield & Murray 1992; Aksnes & 
Franklin 2001; Zhang & Hamilton 2007, 2008), although it could be much larger (Greenberg 
et al. 2008). The stellar value is n °t very important in the case of HD 40307 because 
the planet masses and radii are relatively small (see Eqs. [U-[2]); we assume it is 3 x 10 6 
(Jackson et al. 2008a). 



We must also estimate stellar and planetary radii. We assume that the star's radius 
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follows the empirical relationship found by Gorda Sz Svechnikov (1999). For terrestrial cases 
we scale R p as m^ 27 (Fortney et al. 2007). For gaseous cases, R p is calculated by assuming 
the planet has the same mean density as Neptune: r b = [m b / 'Mjvep) 1 ' 3 R-Nep, where M Nep 
and Rwep are the mass and radius of Neptune, respectively. Note that this assumption is 
consistent with observations of the transiting planet G J 436 b (Deming et al. 2007; Gillon et 
al. 2007; Jackson et al. 2008a). Furthermore, our analysis depends on the age of this system, 
which we estimate as 6.1 ± 1.6 Gyr old (see Appendix). 

3. Results 

Table 1 lists a set of values for the masses and orbits of the planets in HD 40307 
computed by Mayor et al. (2008) with all eccentricities treated as free parameters. P is the 
orbital period and T p is the time of periastron passage. As e values <10 -2 are not currently 
measurable in radial velocity data (the uncertainty is larger than the nominal values), Mayor 
et al. (2008) prefer a solution in which all eccentricities are set to zero (and longitudes of 
periastron vj are therefore undefined). In this case the residuals dropped slightly, but the 
other parameters remained the same (with T p now the time of passage through longitude 
zero). Here, we consider two possibilities for this system, one in which the eccentricities have 
the reported non-zero values listed in Table 1 and one in which the eccentricities are all zero. 

Over timescales much shorter than tidal evolution, interactions among the planets cause 
periodic variations in orbital elements. Fig. [1] shows the variations of the eccentricities, 
produced by our N-body models. In the top panel we model the oscillation of the orbits 
using the reported, non-zero eccentricities. In the bottom panel, we show the evolution 
assuming all eccentricities are initially zero. For these cases we set w b = Wd = and 
w c = 180°. (Different choices of the cc's can result in about a factor of two difference in the 
eccentricity oscillation amplitudes, but this difference does not affect our conclusions.) If we 
assume the eccentricities are initially nonzero, then they remain below 0.1, with planet b's 
eccentricity as large as 0.07. If all the eccentricities are initially zero, then they all quickly 
grow to nonzero values. Even if eccentricities were fully damped by tides, mutual interactions 
would keep them > 0. The minimum perturbed eccentricity is thus ~5x 10 -4 . 

We have modeled the long-term effect of tides on planet b by integrating Eqs. (JT] [2]) back 
in time for various assumed values of its current eccentricity (eo) (Fig. [2]). (The "Alternative 
Model" is explained below.) For these three cases, a& jumps up in values at about the same 
time as e& gets large. So, for example, with eo = 0.008, a& would have been at the location of 
the 2:1 resonance with planet c less than 1 Gyr ago (dotted line in the bottom panel of Fig. 
[2]). If planet b crossed the resonance, both planets' eccentricities would have been pumped 
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up because it is a divergent crossing (Hamilton 1994; Chiang, Fischer & Thommes 2002; 
Zhang & Hamilton 2007, 2008). Such a crossing would likely have destabilized the system, 
or at least prevented the system from appearing as packed as it is today This process would 
have been similar to models of the 2:1 resonance crossing of Jupiter and Saturn in our Solar 
System (Gomes et al. 2005), which significantly spread out the giant planets. For HD 40307, 
such a history is unlikely given the current orbital architecture. Thus, if we assume that 
the resonance crossing could not have happened during this system's history, either (i) tides 
must have damped eo to its minimum possible value several billion years ago, (ii) the system 
must be younger than 1 Gyr, or (iii) Q' b must be larger than 500. 

While using eo = 5 x 10~ 4 as an initial condition for the tidal evolution gives resonance 
crossing 2 Gyr ago, a different history is possible. The current value of e for this planet would 
be ~ 5 x 10 -4 even if the system was fully-damped by tides much earlier. Thus, we have no 
way to constrain how long the system has been in this state; the eo = 5 x 10~ 4 curve in Fig. 
[2] could be shifted by any amount to the right and still be a possible model of the system's 
history. For example, the curve labeled "Alternative Model" has been shifted so that the 
system formed inside the 2:1 resonance 6.1 Gyr ago (its estimated age; see Appendix), then 
was fully damped within 2 Gyr, and remained with the minimum perturbed eccentricities 
ever since (with such small eccentricities, da/dt pa for the intervening 4 Gyr). This scenario 
avoids a past resonance and permits super-Earths (Q' p ~ 500) in the HD 40307 planetary 
system. 

Next we calculate the tidal heating of the planets, still assuming terrestrial planet pa- 
rameters (Q' p = 500, density of 5 g/cm 3 ). For eccentricities of 5 x 10~ 4 , the heating is 
4 x 10 -3 W/m 2 for planet b and 10~ 5 W/m 2 for planet c, assuming radii of 1.07 timeslO 4 
km and 9.4 x 10 3 km, respectively. On the per-unit-surface-area basis (relevant for surface 
characteristics), these values are considerably lower than the Earth's radiogenic heating, 0.08 
W/m 2 . If the planets have the non-zero eccentricities listed in Table 1, the current surface 
heat fluxes for these two planets would be hb = 1 W/m 2 and h c = 0.4 W/m 2 , values that 
lie between the Earth (~ 0.08 W/m 2 ) and Io (~ 2 W/m 2 ). (Note that the Q' value for Io is 
probably similar to the Earth's [Yoder 1995; Aksnes & Franklin 2001]). Even though hb is 
less than the heating on Io, it is probably close enough to result in intense volcanism. 

The short-term periodic variations of the heating due to the periodic variations in ec- 
centricities (from Fig. [1]) are also interesting. The top two curves in Fig. [3] show the changes 
in h that correspond to the orbital changes in the top panel of Fig. [1] h c oscillates between 
0.2 and 0.5 W/m 2 , always more than the Earth's heating. However, hb reaches well over 
10 W/m 2 , and maintains that rate most of the time. Although the heating rates oscillate, 
the periods of oscillation are much shorter than heat transport (predominantly mantle over- 
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turn) in the Earth (see e.g. Davies 1999), so the planet's surface flux would probably be the 
average heating rates (~ 10 W/m 2 ), still several times that of super- volcanic Io. 

Assuming the system has damped to the minimum eccentricities allowed by mutual 
perturbations (Fig. [1] bottom), the heating rates, are much lower; h\, is shown in the lowest 
curve in Fig. [31 Here the heating rate oscillates by many orders of magnitude, but remains 
much less than the Earth's. Planet c's tidal heating in this damped case is always < 10~ 4 
W/m 2 . Therefore, if the eccentricities have damped to minimum values, these planets would 
have heat fluxes much less than the terrestrial planets in our Solar System. 

For the planets orbiting HD 40307 to be rocky, they must have begun tidal evolution 
with low eccentricity orbits and with planet b interior to the 2:1 resonance with b. How- 
ever, it is unlikely that the planets acquired most of their mass in such a configuration, 
as that scenario predicts an implausibly large pre-planetary nebula. According to Kuch- 
ner (2004), such in situ formation requires a primordial disk with surface density profile 
E(r) = 6379(r/0.1AU)" - 925 gem" 2 . This disk contains 21.6 M e inside 0.15 AU, probably 
15-100 times more than our Solar System had (Weidenschilling 1977; Hayashi 1981) @ This 
result is consistent with the assumption that HD 40307's gas-to-dust ratio is half the solar 
nebula's because [Fe/H] = -0.31, in which case its disk would have been 30-200 times more 
massive than the solar nebula. These values represent a disk mass comparable to the stellar 
mass, and far in excess of the typical star-disk mass ratio of 1% (Andrews & Williams 2005). 
We conclude that these planets did not form in situ; they must have formed further out and 
migrated in prior to the dispersal of the gas disk. 

So far we have considered the implications of rocky planets. Could they instead be 
mini-Neptunes with a thick gaseous envelopes and Q' = 10 5 ? If the current e& has the value 
given in Table 1, the tidal history (shown by dashed lines in Fig. [2]) would not have included 
dangerous resonance crossings in the last 6.1 Gyr. Furthermore, such planets would have 
heating fluxes several orders of magnitude smaller than terrestrial planets, and their internal 
structures would probably not be significantly affected by tidal heat. If, however, the mass 
(and hence radius) of planet b is significantly greater than the observational minimum, then 
there may have been more evolution. We find that if = 15 M e , then HD 40307 b would 
have been at the 2:1 resonance 6.1 Gyr ago. Such a mass corresponds to an orbit inclined by 
75° to the line of sight. In other words, planet b's orbit must be more than 15° from face-on. 

If the planets are gaseous, did they form in situ, or did they migrate in from further out? 
If the planets formed in situ, the pre-planetary nebula would have required about half as much 



2 We assume that E = E (r/1AU)^ :E , where (E Q ,x) = (7.75 gcm~ 2 , 1.5) or (5.895 gem' 2 , 1.0) and extend 
the disks in to r = 0. These disks are calibrated to the MMSN and contain 5 M ffi from 0.4-4 AU. 



- 8- 



mass as for the rocky-planet case described above (assuming the cores' masses are roughly 
equal to the envelopes' masses), but that value is still improbably large. More likely, as with 
rocky planets, the planets could have formed further out and migrated in via interactions 
with the disk (Lin et al. 1996). Theoretical models of this phenomenon (e.g. Snellgrove et 
al. 2001; Lee & Peale 2002) suggest resonance capture could be an outcome of this process, 
but resonances are not observed (Mayor et al. 2008). However, such a commensurability 
could have been destroyed by subsequent mergers, scattering or turbulence (Terquem & 
Papaloizou 2007; Adams, Laughlin & Bloch 2008; Lecounet et al. 2008). Therefore the 
migration scenario is consistent with the observed orbits. We conclude that if the planets 
are mini-Neptunes, they likely formed at larger radii and migrated in via disk torques to 
their present orbits. 

The terrestrial and gaseous planet models considered above are not a complete explo- 
ration of parameter space. For example, the responses of the planets to tides are encapsulated 
in Q', a notoriously uncertain parameter, even in the Solar System. To address this uncer- 
tainty, we solved Eqs. (P - |2J) for a range of values of Q' and e (1500 < Q' b < 3500 and 
10~ 4 < eo < 0.01) to determine how long ago the resonance crossing would have occurred. 
If it occurred more than 6.1 Gyr ago, the case is allowed, but if the crossing occurred within 
the last 6.1 Gyr, then the case is forbidden. For ej, < 10~ 3 , these restrictions do not strictly 
apply since the eccentricity may have damped to those values an arbitrarily long time ago. 
If, however, e& is found to be greater than 10 -3 , then it would constrain Q' b to be >2200. 



4. Conclusions 

The small minimum masses of the planets orbiting HD 40307 are tantalizingly close 
to masses of the rocky, "terrestrial" planets in our Solar System. By considering the dy- 
namical features and history of the system, we have determined implications of their being 
predominantly rocky or gaseous. We find both possibilities are consistent with the observa- 
tions, but the likelihood that they are terrestrial depends on the actual values of the current 
eccentricities, which are at or below the detection threshold. 

Mayor et al. (2008) report two sets of e values, either all ~ 0.01 or all zero. The latter 
case is ruled out by mutual perturbations between the planets (see Fig. [[]). Instead, the 
lowest possible eccentricities are ~ 10 -3 . If the values are that small, tides probably damped 
them from higher values, and they may have reached this state recently, or at any time in the 
past. These ambiguities preclude a definitive assessment of the composition of the planets. 
However, if or when the eccentricities are measured more precisely, such a determination will 
be possible through tidal analysis. 
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If the planets are terrestrial (i.e. Q' ~ 500), then either (a) the system must be less 
than 2 Gyr old, which is unlikely (see Appendix), or (b) tidal evolution began with planet 
b just inside the 2:1 resonance with c, and with modest eccentricity (< 0.3). In either case, 
the planets must have formed at larger distances and migrated inward. In case (a), the 
inner planet is a "super-Io", with intense volcanism, a type of rocky body suggested by 
Jackson et al. (2008c). In case (b), the planets could be terrestrial- like bodies: rocky, with 
thin atmospheres and modest volcanism. 

If the planets are "mini-Neptunes" , then the orbital history, formation scenarios, and 
internal structures are all consistent with previous models of such bodies. The only constraint 
tidal evolution can provide is that i > 15° (assuming Q' b = 10 5 ). As this constraint is less 
stringent than the requirements for a terrestrial body, our analysis favors the mini-Neptune 
model somewhat. 

In order to distinguish the super-Earth and mini-Neptunes through dynamical analyses, 
the eccentricities need to be determined to within at least 10~ 3 . Such high precision is only 
measurable for planets that undergo secondary transits, but not even a primary transit has 
(so far) been reported for this system. Photometric observations of the transit would allow 
a calculation of the planet's radius, but this value is not enough to determine if the plan- 
ets are super-Earths because of composition-radius degeneracies exist in current models of 
small-mass exoplanets (Adams, Seager & Elkins-Tanton 2008; Raymond, Barnes & Mandell 
2008). Therefore dynamical models of the system's history may be the most effective way 
to determine the composition of the planets. 

We also note that if planet b is gaseous, then it may be undergoing evaporative mass 
loss (Baraffe et al. 2004; Hubbard et al. 2007). Indeed, Raymond et al. (2008) showed that a 
25 M ffi planet at 0.05 AU could be evaporated to its core in 4-5 Gyr, although this depends 
on the star's X-ray history. If planet b had its atmosphere removed in the past, then the 
tidal models presented in § 3 may be inadequate as we have not considered time- varying Q' b . 
Furthermore, even if HD 40307 b is gaseous today, then it may one day be reduced to a solid 
core, making it terrestrial. 

Our analysis also admits the possibility that HD 40307 b is an exotic planet, unlike 
any in our Solar System, with end-member properties ranging from a volcanic super-Io to 
mini-Neptunes. Q' values may range from ~ 10 3 to ~ 10 4 , perhaps resulting from an unusual 
internal structure. 

The dynamical properties of planetary systems may be used to constrain the physical 
properties of exoplanets as Fig. H] demonstrates. As more potentially terrestrial-like planets 
are detected, dynamical analyses will continue to play a role in constraining their physical 



and orbital properties. 
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Table 1. Masses and orbits for the HP 40307 planets 



Planet m ( M©) P (d) a (AU) e T p ( JD) 



b 


4.2 


4.3115 


0.047 


0.008 ±0.065 


2454562.77 


c 


6.8 


9.62 


0.081 


0.033 ±0.052 


2454551.53 


d 


9.2 


20.46 


0.134 


0.037 ±0.052 


2454532.42 



a Set to zero by Mayor et al. (2008) 
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Fig. 1. — Top: Eccentricity evolution of HD 40307 b, c, and d assuming nominal parameters from 
Table 1. Bottom: Eccentricity evolution assuming all e values are initially zero (the timescale here 
is shorter than above because the oscillation period is much shorter and was chosen to illustrate 
the maximum values of e^). 
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Fig. 2. — Top: History of eb for different orbital and physical assumptions. Solid curves assume 
Q' b = 500 (rocky super-Earths), dashed curves assume Q' b = 10 5 (gaseous mini-Neptunes). The 
assumed value of the initial eccentricity of planet b, eo, is indicated on each curve. The "Alternative 
Model" curve assumes that planet b formed with e w 0.25 and a& just interior to the 2:1 resonance 
with planet c. The value of in this case damped to its minimum perturbed value and remained 
there. Bottom: Evolution of a b (solid and dashed lines) for the same cases and in the same 
sequence as above. For reference the location of the inner 2:1 resonance with planet c (labeled "2:1 
Resonance") is shown by the dotted line. Note that we assume that any history prior to the 2:1 
crossing is unphysical. 
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Fig. 3. — Tidal heating fluxes of terrestrial- like planets in HD 40307. The upper two curves 
correspond to the changes in e shown in Fig. Q] (top), which used initial conditions from the nominal 
orbital elements in Table 1. The curve labeled "b-damped" corresponds to the changes in e shown 
in Fig. [I] (bottom), in which orbits are initially circular. 
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Fig. 4. — Values of Q' b and that predict resonance crossing less than 6.1 Gyr ago (forbidden 
region) and those that don't (allowed region). This plot demonstrates how the orbital properties 
of the system can constrain the physical properties of planet b. Note that the Earth's Q' value is 
~ 500 and Neptune's is ~ 10 5 . 
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Appendix: The Age of HD 40307 

Based on two techniques that use the rotation rate and the magnetic activity strength 
respectively, we estimated the age of HD 40307. Dwarf stars with spectral types ranging 
from late-F to mid-M have both a radiative and a convective zone, the interface of which 
is thought to be responsible for magnetic field generation and the heating of the upper 
atmosphere (Parker 1993; Ossendrijver 2003; Thompson et al. 2003), a phenomenon that 
gives rise to magnetic activity. These stars begin their lives rotating quickly but slowly 
lose angular momentum over time via solar-type winds. This loss of angular momentum as a 
function of stellar age has been observed and can be quantified (e.g. Skumanich 1972). Using 
the age-rotation relation for a star with a radiative-convective zone interface, we calculated 
an age of 6.7±2.0 Gyr based on a 48 day period, and using the 30% quoted uncertainty from 
Barnes (2003; I sequence). 

Since the rotation rate is linked to magnetic field generation and the subsequent chromo- 
spheric activity in solar-type stars, a reduced rotation rate will result in less magnetic activity. 
Using the cluster derived age-activity relations from Soderblom Duncan & Johnson (1991), 
we calculate an age (averaging the 3 relations) of 5.4±1.6 Gyr, based on an R' HK value of 
-4.99 (Noyes et al. 1984). These results are confirmed using the age-activity relation of 
Donahue (1998), which also yield an age of 5.4 Gyr. Recently, Mamajek & Hillenbrand 
(2008) refined the age-activity relations for F-K dwarfs using additional stellar clusters with 
improved age estimates. Their relations yield an age of 6.3 ± 0.9 Gyr for HD 40307. The 
Donahue (1998) study uses coeval binary systems to quantify the uncertainty in these rela- 
tions. For the early K-type dwarfs, the age-discrepancy in binary pairs is smaller than the 
uncertainty in the age-activity relation, confirming an uncertainty of ~ 1.6 Gyr. 

Combining the results from the rotation and activity analysis, we estimate the age of HD 
40307 to be 6.1±1.6 Gyr. 



